Subwavelength position sensing using nonlinear feedback and wave chaos 
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We demonstrate a position-sensing technique that relies on the inherent sensitivity of chaos, where 
we illuminate a subwavelength object with a complex structured radio-frequency field generated 
using wave chaos and a nonlinear feedback loop. We operate the system in a quasi-periodic state 
and analyze changes in the frequency content of the scalar voltage signal in the feedback loop. 
This allows us to extract the object's position with a one-dimensional resolution of ~ A/10,000 and 
a two-dimensional resolution of ~ A/300, where A is the shortest wavelength of the illuminating 
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PACS numbers: 05.45.Gg, 05.45.Mt, 84.30.Ng 

Diffraction, a property of electromagnetic (EM) waves, 
blurs spatial information less than the wavelength A of 
an illuminating source and hence limits the resolution of 
images. Over the past decade, techniques have been de- 
veloped that overcome this diffraction limit using super- 
lenses made from negative- index media [l|, Q, super- 
oscillations 0], and nano-structures with surface plas- 
mons [3, HI]. Other methods use fluorescent molecules 
that serve as subwavelength point markers where 
imaging is enabled by sensing the position of the markers. 

In this Letter, we describe a new super-resolution tech- 
nique that senses the position of an object by combin- 
ing two concepts: nonlinear delayed feedback and wave 
chaos. The system uses radio frequency (RF) EM waves 
in a closed feedback-loop through a wave-chaotic cavity. 
Self oscillation in the feedback occurs when the loop gain 
exceeds the loop losses; no RF field is supplied by an 
external source. We include a nonlinear element (NLE) 
in the feedback loop to create a system with complex 
(non-periodic) dynamics. The resulting EM oscillations 
provide the illumination source for our position sensor. 

Our work extends the Larsen effect, a known phe- 
nomenon where positive audio-feedback between a mi- 
crophone and audio amplifier results in periodic acoustic 
oscillations. The frequency of oscillation, known as the 
Larsen frequency, is highly dependent on the propaga- 
tion paths of the acoustic wave. A perturbation to these 
propagation paths shifts the Larsen frequency fioj . For 
our super-resolution position-sensing system, we exploit 
the sensitivity of quasi-periodic EM frequencies. 

In our experimental system, the NLE is an in- 
put/output circuit based on the design from Ref. [llj ]. 
We use an aluminum two-dimensional (2-D) quarter- 
stadium-shaped RF cavity for a wave-chaotic scattering 
scene [l2|. As shown in Fig.QJ,, the EM field emanating 
from the cavity is fed into a nonlinear circuit through a 
broadband (20 MHz - 2 GHz) receiving antenna (RX). 
The output of the circuit is fed back into the cavity 
through an identical transmitting antenna (TX), creating 
a closed feedback loop. Inside the cavity is a subwave- 
length dielectric object. 
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FIG. 1: (a) Experimental setup. Nonlinear circuit consist- 
ing of a transistor-based (BFG620) NLE, low-noise ampli- 
fiers (LNA, Mini-Circuits ZX60-4016E and Picosecond Pulse 
Labs 5828-108), and a low-pass filter. A bias voltage Vbias 
tunes the nonlinearity. The output voltage V ou t is measured 
with a 8-GHz-analog-bandwidth 40-GS/s oscilloscope (Agi- 
lent DSO80804B). Feedback passes through a cavity with a 
dielectric object (2 cm x 4 cm water-filled container) that is 
positioned in 2-D using Thorlabs (LTS150) and Zaber Tech- 
nologies (TLSR150B) translations stages, (b) Cavity pulse- 
response. We inject a 0.1 ns EM pulse of amplitude 1.5 V 
through the TX antenna and measure at the RX antenna. 
The time at which the radiation arrives at RX is a measure 
of the path length of EM energy through the cavity. 



The complex field inside of the cavity interacts mul- 
tiple times with the object due to reflections from the 
cavity's walls (illustrated by a complex ray path in Fig. 
UK) and from many passes of the RF signal through the 
nonlinear feedback loop. Due to these multiple interac- 
tions, small object movements change the structure of 
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the field. These changes alter the dynamical state of the 
system. The output of the nonlinear circuit is filtered 
such that its maximum frequency is 2 GHz, and thus the 
RF signal has A > 15 cm. 

The NLE in the EM feedback loop induces quasi- 
periodic oscillations with multiple incommensurate fre- 
quencies in the output voltage V ou t of the nonlinear cir- 
cuit. As the object moves inside the cavity, the frequen- 
cies of the quasi-periodic oscillations shift independently 
and provide a unique fingerprint of the object's location 
in 2-D. Thus, we map the position of the object in both 
the x and y directions by monitoring changes of a single 
scalar voltage V ou t- 

Before describing our results, we first characterize our 
wave-chaotic cavity using its pulse response. Shown in 
Fig.QJ), our cavity produces a complicated pulse response 
(typical of wave-chaotic systems) with a quality factor 
Q = 174 at a frequency of 1.77 GHz (the most promi- 
nent frequency in the quasi-periodic oscillations). As a 
result, broadcasting a continuous-wave signal into this 
cavity forms a complex interference pattern for each con- 
tained frequency (generic cavities tend to display such 
wave chaos; only cavities with a high degree of symme- 
try display simple interference patterns) [l 31 ] . 

The pulse response of a wave-chaotic environment has 
been exploited to sense the appearance of an object in a 
scattering medium jl4| or the location of a perturbation 



on the surface of a scattering medium [15j. These tech- 
niques rely on measuring changes to the pulse response 
and have demonstrated a spatial sensitivity of ~ \. Our 
own work is inspired by these achievements, where we 
use a continuous-time nonlinear feedback loop to achieve 
deep subwavelength position resolution. 

Conventional oscillators using time-delayed nonlinear 
feedback use a nonlinear element whose output is ampli- 
fied and coupled back to the input through a single feed- 
back loop that delays the signal by a fixed amount. These 
systems can display a variety of behaviors including peri- 
odic oscillations, quasi-periodicity, and chaos. Oscillators 
using time-delayed feedback have been designed using 
high-speed commercial electronics or lasers to generate 
complex signals with freque ncy bandwidths that stretch 
across several gigahertz [16h18| . 

Thus, the system shown in Fig.QJi combines the sensi- 
tivity of a dynamical state from a high-speed nonlinear- 
feedback oscillator with the sensitivity of the EM field in 
a wave-chaotic cavity. The time delays of the feedback 
in this system are the propagation times for the EM en- 
ergy to transmit through the cavity, rather than a single 
time-delay. The values of the delays and their respective 
gains form a continuous delay distribution (proportional 
to the cavity's pulse response) that is uniquely defined for 
each position of the enclosed object. Due to the nonlin- 
ear feedback, the system's dynamics are highly sensitive 
to changes in this distribution of delays. Measuring the 
scalar variable V ou t , we monitor dynamical changes in the 
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FIG. 2: Bifurcation Diagram. We store and analyze time 
series of V ou t at each object position along a 1-D path in the 
x direction. The local maxima of each time series is plotted 
as a function of object position x. 



system and sense the object's movements. 

We first demonstrate this idea qualitatively along a 
one-dimensional (1-D) object path. We fix Vbias in the 
NLE to exhibit periodicity at x = mm and measure the 
time evolution in V ou t for object positions x = mm - 
12 mm in 10 /im steps. The system changes between pe- 
riodicity (P) from x = mm - 1.4 mm, quasi-periodicity 
(QP) from x = 1.4 mm - 8 mm, and two different time- 
evolving chaotic states (CI and C2) from x — 8 mm - 
9.8 mm and x — 9.8 mm - 12 mm, respectively Chaotic 
state CI contains chaotic- like breathers and C2 exhibits 
a relatively flat bandwidth from 20 MHz - 2 GHz. 

The observed dynamical changes fall into one of two 
categories: an abrupt change in the dynamical state 
(known as a bifurcation) or small shifts in the frequency 
components and amplitudes of V ou t- A bifurcation di- 
agram illustrates the qualitative dynamical changes in 
Fig. [2j Our results show dynamical changes from sub- 
wavelength movements of a subwavelength object. 

To go beyond the qualitative detection of movement, 
we tune Vbias so that V ou t is in a quasi-periodic state 
(QP in Fig. [2]) for all object positions of interest. The 
incommensurate frequencies of a QP state are not phase- 
locked and hence can shift independently with respect 
to object translations. In addition, incommensurate fre- 
quencies help eliminate interference nodes (blind spots) 
of the illuminating EM waves in the cavity, where each 
frequency has a complex interference pattern that covers 
the blind spots of another. An example time series and 
frequency spectra for a fixed object position are seen in 
Fig.EK and Fig. [3j}, respectively. Recall that, though the 
system's dynamics are QP in time, the EM energy inside 
of the cavity is chaotic. 

Tracking the object entails measuring shifts in the QP 
frequency components. In Fig. [3b, we highlight two peaks 
in the spectrum at frequencies denoted by f\ and f^. 
The frequency harmonics at (/2 — fi) and (/2 + /i) are 
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FIG. 3: (a) Temporal evolution of a typical quasi-periodic 
state used for position sensing, (b) Logarithmic power spec- 
tral density (PSD) of this signal. 

used to improve the signal-to-noise ratio (SNR) of these 
frequencies. Averaging independent measures of /i and 
f 2 reduces statistical errors and increases their SNR. To 
follow changes in the frequencies of V ou t with high preci- 
sion, we use a nonlinear least-squares-fit to a model for a 
four-tone QP signal 0, resulting in a 2.4 kHz frequency 
resolution (approximately 0.5% of the total observed ex- 
perimental frequency shifts). 

To demonstrate 1-D position sensing, we translate the 
object along the path x = mm - 5 mm while y — 2.5 
mm. We then translate the object along an orthogonal 
path y — - 5 mm while x = 2.5 mm. Shown in Fig. |4^, 
and Fig. HJj, the measured frequency shifts A/i and Af 2 
are plotted for 1-D paths along the orthogonal x and y 
directions, respectively. We separately fit A/i and A/ 2 
in the x and y directions with second order polynomials 

aiA/i(x) + a 2 Af 2 (x) = c + c x x + c 2 x 2 , (1) 

hAhiy) + b 2 Af 2 {y) = d Q + d lV + d 2 y 2 . (2) 

We optimize the coefficients aj and Cj (6j and di) using 
a nonlinear least-squarcs-fit to a model for the object 
position. The root-mean-square (RMS) errors for the 
frequency shift map is 1.45 kHz (0.86 kHz) along x (y). 
By inverting these maps, we calculate the measured ob- 
ject positions. The RMS error between the actual and 
measured positions is 9.2 ^m (23.7 /im) for x (y), which 
demonstrates a resolution of ~ A/10, 000 along orthogo- 
nal 1-D directions (recall A > 15 cm). 

Tracking the object's position in both the x and 
y directions simultaneously requires two independently 
changing observables. In our system, we observe a single 
scalar variable V^ ut that oscillates with primary frequen- 
cies /1 and f 2 . We fit the frequency shifts Af±(x,y) and 




FIG. 4: Frequency shifts A/i and A/2 of the QP state as 
the object translates along the (a) x direction, where the fit 
from Eq. ([TJ yields a\ = 2.8 mm/kHz, a 2 = -8.7 mm/kHz, 
Co = 27,113.0 mm, ci = -439.2 and C2 = 1.5 mm -1 , and (b) 
along the y direction, where the fit from Eq. ([2]) yields b 1 = 
0.8 mm/kHz, b 2 = -0.5 mm/kHz, d = -73,452.0 mm, di = 
242.0, and d 2 = -2.0 mm -1 , (c) Frequency shifts Aft(x,y) 
and Af2(x,y) for object translations (x,y) in a 5 mm x 5 
mm grid of positions. The planar fits from Eqs. (0 and @ 
yield ai = -84.68 kHz/mm, q 2 = -15.20 kHz/mm, ft = -56.74 
kHz/mm, ft = -14.75 kHz/mm, ei = 11.72 MHz and e 2 = 
2.43 MHz. The measured determinant |aift — 0:2ft | = 386 
kHz 2 /mm 2 with an error of 5.8 kHz 2 /mm 2 . 

A/a(a;, y) for object positions (x,y) in a 5 mm x 5 mm 
area (Fig. |4j^) and approximate them as planes 

&h(x,y) = a 1 x + firy + ei, (3) 

A/ 2 (x, y) = a 2 x + (3 2 y + e 2 . (4) 

Using Cramer's rule, we show that \ai(3 2 — a 2 (3i\ ^ to 
verify the planes are linearly independent in this area and 
allow us to simultaneously measure x and y coordinates. 

In the 1-D case, we have the freedom to optimize the 
fitting parameters in Eqs. |T]) and ^ for x and y sepa- 
rately. In the 2-D case, all of the fitting parameters o^, 
fti, and e, in Eq. and Eq. (01 are present in the so- 
lutions for both x and y. Thus, we cannot optimize the 
fits in the x and y directions separately and instead use 
the fitted planes for our frequency maps. 

This constraint, combined with the approximation that 
these surfaces are planar, limits our 2-D resolution. A 
planar fit of A/i(x, y) (Af 2 (x, y)) gives a RMS frequency 
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error of 4.17 kHz (7.26 kHz) and a RMS position error of 
370 fj,m (650 /xm) for x (y), yielding a 2-D resolution of 
~ A/300. Higher order fits do not improve the resolution 
due to noise in our measurements. This 2-D frequency- 
mapping serves as the calibration for objects of this shape 
and must be reacquired for different shaped objects. 

For comparison, a scanning near-field microwave mi- 
croscope uses RF frequency shifts to achieve subwave- 
leng th sensitivity (~ A/750,000) of near planar surfaces 
[20j. In contrast, our system uses nonlinear feedback 
to internally generate multiple independent frequencies 
and measures multiple degrees-of-frcedom using a single 
scalar variable. Moreover, it uses a stationary pair of an- 
tennas to extract 2-D spatial information of a 3-D object, 
making it free of mechanically-moving parts. 

We conjecture that our method can be implemented 
using EM waves in the visible part of the spectrum. 
Semiconductor lasers with time-delayed optical feedback 
are known to display complex dynamical behaviors in 
which the output intensity varies in time, including quasi- 
periodicity [2ll423j . Furthermore, optical wave chaos has 
been demonstrated using optical cavities (24l - [26| . We en- 
vision a completely optical version of our technique where 
a laser receives feedback from a wave-chaotic optical cav- 
ity. Such a system will be capable of tracking an object 
on a sub-nanometer scale. 

Understanding the full potential of this method will 
require studies in both wave chaos and nonlinear dynam- 
ics. Our results suggest that one can position sense in 3- 
D using a QP state with three independent frequencies. 
This type of dynamical state is possible in our system 
but requires further study to create a QP state for a 3-D 
volume of interest. More independent observables could 
also be introduced into the system using two or more 
feedback loops external to the cavity, where each loop is 
independently band-limited to prevent cross talk. 

In the future, we see several options to improve the sys- 
tem's resolution. Increasing the number of frequency har- 
monics through nonlinear mixing gives additional mea- 
sures of the independent modes and improves the sys- 
tem's SNR. Also, the cavity Q is proportional to the 
number of interactions between the subwavelength ob- 
ject and the EM energy inside of the cavity, and thus the 
resolution of this technique should also scale with Q. 

We believe that our system will have applications be- 
yond position sensing. Subwavelength scatterers are of- 
ten treated as point-like objects; our approach is sensitive 
to the shape and orientation of the subwavelength scat- 
tcrcr. Also, similar to [lol ]. analyzing dynamical states 
can monitor changes in the EM properties of materials 
in the cavity. 

To the best of our knowledge, our approach is the first 
to measure multiple spatial degrees-of-freedom on a sub- 
wavelength scale using a single scalar signal. Using a 
QP analog of the Larsen effect, we combine a nonlin- 
ear feedback oscillator with multiple EM reflections in a 



scattering environment to exploit the inherent sensitivity 
of wave chaos, adding an alternative to the short list of 
super-resolution techniques. 

We gratefully acknowledge Zheng Gao with help in de- 
signing NLE and the financial support of the U.S. Office 
of Naval Research grant # N000014-07-0734. 
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